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Abstract 

We develop a new semi-dynamical method to study shock revival by neutrino heating in core- 
collapse supernovae. Our new approach is an extension of the steady state approach with the light 
bulb approximation, which has been widely used in the supernova community for the investigation 
of the criteria for successful supernova explosions. In this approach we get rid of the steady- state 
condition and take into account deviations from it instead. We have in mind the scenario that not the 
critical luminosity but the critical fluctuations generated by hydrodynamical instabilities such as SASI 
and neutrino- driven convection in the post-shock region determines the onset of shock revival. We 
validate our new approach by some test calculations, confirming that it indeed captures the essential 
evolution of shock wave qualitatively. We then apply the method in various conditions to find that 
there is a critical fluctuation for shock revival. The critical fluctuation can be well fit by the following 
formula: f cr i t ~ 0.8 x (Mi n /lAMo) x {1 — (r s / l /10 8 cm)}, in which / denotes the normalized amplitude 
of pressure perturbation. We also find that the critical fluctuation decreases with the shock radius, 
whereas it increases with the mass of the central core. We discuss the possible implications of our 
results for 3D effects on shock revival, which is currently controversial in the supernova community. 
Subject headings: supernovae: general — neutrinos — hydrodynamics 



1. INTRODUCTION 

Decades of research on core-collapse supernovae (CC- 
SNe) have not accomplished full understanding of the ex- 
plosion mechanism yet. The fundamental question in su- 
pernova theory is how a shock wave propagates through 
the stellar core, mantle and envelope and finally dis- 
rupt a star. It is widely recognized that the prompt 
shock wave, which is generated by core bounce, expe- 
riences stagnation at r ^ 100km owing to the energy 
losses by photodissociations of heavy elements as well as 
neutrino cooling and ram pressure of accreting matter. 
For successful explosions, the stalled shock wave needs 
re-invigoration one way or another. The most promis- 
ing energy source is neutrinos diffusing out of a proto- 
neutron star. Only a percent of the total amount of 
energy neutrinos are carrying is sufficient to given the 
canonical explosion energy of supernova. Although it is 
still unclear whether the energy transfer from neutrinos 
to eject a is indeed large enough, the neutrino heating 
is c urrently the most fav ored mechanism of CCSNe (see 
e.g. iKotake et al.l (|2012l ) for a recent review). 

In the the last five years, we have seen some ever more 
evolving multi-dimensional num erical simulations pro - 



It was iBurrows fc Goshvl (|l993l ) who took the lead in 
such an approach. They examined a sequence of steady 
accretion flows through a standing shock wave onto a 
proto-neutron star to discuss the criterion for shock re- 
vival. They varied the electron- type neutrino luminos- 
ity (L Ue ) and mass accretion rate (M) as free parame- 
ters and revealed that there is a critical luminosity for 
a given mass accretion rate, which was defined to be 
the threshold, above which no steady shocked accretion 
flow obtains. They argued that the critical luminos- 
ity marks the trigger point for shock revival. This ap- 
pr oach was later extended to rotational configurations 
bv lYamasaki fc Yamad a (2005) and linear stability anal- 
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vsis was also applied ([Yamasaki fc Yamadal 120071 ) . The 
latter authors pointed out that the critical luminosity 
could be smaller if one takes into account hydrodynam- 
ical instabilities such as radial overstabilization modes, 
whic h were also observed in dynamical simulations (see 
e.g. lOhnishi et al.l ([2006)) as well as non-radial modes. 
Although the approach certainly has limitations in repro- 
ducing full complexity of explosion mechanism, the crit- 
ical luminosity has become one of the most useful mea- 
sures for shock revival and some ne w analyses have been 
publis h ed in recent years (see e .g. iPejcha fc Thompson! 
(|2012f ): lKeshet fc B albersl (|2 01 2f ) and references therein). 

The existence of the critical luminosity has been also 
demonstrated in multi-dimensional numerical simula- 
tions with simplified treatments of neutrin o transfer such 
as th e gray or light bulb approximation (lOhnishi et al.l 
120061: Ew akami et al. 2008; Mur phv fc Burrows! 120081 : 
iNordhaus et al.ll2010l : lHanke et al.ll201lf ). There is a wide 
consensus at present that multi-dimensional effects such 
as th e standing acc retion s hock instability (SASI ) (see 
e.g. iBlondin et all (|2003[ ): iFodizzo et ail (120071) ) and 
neutrino- driven convection (see .e.g. iMurphv fc Meakinl 



duce expl o sion m o dels (see e.g. [ Burrows et al.l (120061): 



lOtt et all (120081): iMarek fc Jankal (120091): ISuwa et all 



iTakiwaki et al.l (j2012l ): iBruenn et al.l (j2010f ): 
iMueller et al.l (|2012f ). These results still have some dis- 
crepancies among different groups, though, and how 3D 
dynamics affects the outcome is remaining to be inves- 
tigated in detail. Although large-scale numerical simu- 
lations are making a rapid progress toward sufficient re- 
ality, we believe that phenomenological approaches such 
as proposed in this paper can still play an important and 
complimentary role for extracting key physical elements 
from the complex and non-linear dynamics in CCSNe. 
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([201 lh : iMurphv et "all (|2Q12[ ) and references therein) re- 
duce the critical luminosity. This is mainly because the 
advection time scale tends to be longer owing to tur- 
bulent motions, which then leads to longer heating and 
provides a favorable condition for shock revival. In addi- 
tion, the instabilities push the shock wave outward and 
expands the gain region. Hence the exploration of multi- 
dimensional neutrino-heating mechanism is currently the 
main topic of main-stream research on CCSNe. 

In pa rticular, 3D dynamics w ill be one of the central 
issues. iNordhaus et al.l (|2Q1Q[ ) found in their 2D and 
3D experimental simulations that the critical luminos- 
ity is monotonically reduced with increasing spatial di- 
mensions, i.e., 3D dynamics provides the conditions that 
are most fav o rable for shock revival. On the other hand, 
lHanke et al.l (|201ll ) obtained in their similar computa- 
tio ns the results that are at odds with those presented 
by INordhaus et al.l (|2Q1Q[ ) : they found no significant dif- 
ference in the critical luminosity between their 2D and 
3D models. The reason of the difference is not clear at 
the moment, since they made different simplifying as- 
sumptions and employed different numerical techniques. 

In this paper, we address the condition for shock re- 
vival again by a simplified phenomenological approach. 
We do not discuss the critical luminosity, however. In- 
stead we propose to introduce the third parameter, fluc- 
tuations, and discuss the critical fluctuation for shock 
revival. In so doing we make an attempt to construct a 
semi-dynamical model to approximately describe shock 
motions that are induced by fluctuations. This is an 
extension of the commonly-employed steady-state ap- 
proach. One of the drawbacks in the latter approach 
is that by definition it cannot handle the temporal evo- 
lution of shock wave and, consequently, cannot address 
the fate of revived shock wave if it indeed revives at the 
critical luminosity in the first place. The semi-dynamical 
approach can remove these drawbacks and demonstrate 
that there is a threshold of fluctuation amplitudes for a 
given combination of neutrino luminosity and mass accre- 
tion rate, beyond which shock revival occurs and a con- 
tinuous outward propagation of the revived shock wave 
follows, which we think is the new criterion for successful 
explosions. 

The paper is organized as follows. In Section 2, we 
describe the semi-dynamical model. We conduct several 
calculations for validation and then apply the model sys- 
tematically to various conditions in the shock-stagnation 
phase in Section 3. Based on these results we give the 
critical fluctuation amplitudes for shock revival. Finally, 
we discuss the possible implications of the new criterion 
for the multi-dimensional neutrino heating mechanism in 
supernova theory and give conclusions in Section 4. 

2. SEMI-DYNAMICAL METHOD 

We are interested in the phase at several hundred mil- 
liseconds after bounce, in which the prompt shock wave 
is stalled and matter flows through the almost steady 
shock and accretes onto a proto-neutron star. Neutri- 
nos are emitted from the neutrino sphere and heat up 
the gain region. The location of quasi-steady stagnated 
shock wave is determined by the neutrino heating and the 
ram pressure of accreting matter. The onset of shock re- 
vival is determined not by the non-existence of steady ac- 
cretion flows but by some sort of hydrodynamical insta- 



bilities, e.g., radial over-stabilizing o scillations in spheri- 
cal symmetry an d SASI in multi-D ([Ohnishi et al.l 120061 : 
lFernandezll20ll ). In this paper we take a standpoint that 
the onset of shock revival is determined by the neutrino 
luminosity and mass accretion rate as well as fluctua- 
tions by the instabilities. Our semi- dynamical approach 
deals with the transition from the quasi-steady accretion 
to the re-expansion of shock wave as well as the ensuing 
outward propagations of shock wave in ID. Before go- 
ing into detail, we first describe the essence of this new 
approach. 

The semi-dynamical model begins with an addition 
of perturbation to the steady accretion shock by hand. 
Then the shock wave starts to move. We follow the subse- 
quent shock motions not by hydrodynamical simulations 
but by the integration of a simplified equation of mo- 
tion, which is base on the local Riemann problem. By 
local we mean that our formulation considers only the 
neighborhood of shock wave, which enable us to avoid 
simulations. This local approximation is found to be in 
good agreement with the results of full dynamical sim- 
ulations near the shock wave (see Section I2.1.2| 12.21 and 
Appendix [Bj . 

The model calculations are computationally very 
cheap, by virtue of which we can investigate long-term 
(> Is) evolutions of shock wave for a large number of 
models with different backgrounds so that we could ob- 
tain the critical fluctuation amplitudes for shock revival 
very efficiently. Another benefit for the semi-dynamical 
approach is that it can study finite fluctuation ampli- 
tudes, which is in sharp contrast to linear stability anal- 
ysis, in which the stability is considered only for infinites- 
imal perturbations. We stress that even if the shock wave 
is linearly stable, large enough perturbations may trigger 
shock revival, which we will see is really the case in the 
later section. 

2.1. Details of the semi- dynamical model 

The actual calculations of shock motions in our semi- 
dynamical approach consist of a couple of steps, which 
we will describe more in detail in the following section 
(see also Figure [1]). 

2.1.1. Step 1: preparation of steady shocked accretion flows 

The first step is a preparation of initial conditions, 
which are assumed to be steady and spherically sym- 
metric. The procedu re to obtain s uch fl ows is ex- 
actly the same as i n lOhnishi et al.l (|2006h ( see also 
lYamasaki fc Yamadal (|2005f )) and the details are given 
in Appendix lAl 

2.1.2. Step 2: addition of perturbations 

After setting up the steady shocked accretion flows, we 
put some radial perturbations, which are not necessar- 
ily small, to the initial conditions. Although the actual 
form of perturbation is rather arbitrary, we choose in this 
study to give a finite velocity to the standing accretion 
shock wave by adjusting post-shock quantities so that 
the Rankine-Hugoniot relation would be satisfied for the 
given shock velocity, v s h(to)- In this expression, to de- 
notes the initial time. 

2.1.3. Step 3: displacement of shock wave 
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Fig. 1. — Schematic pictures of individual steps in the semi-dynamical model. The horizontal axis in each panel denotes the radius 
whereas the vertical axis represents the density. See the text for more detailed explanations of each step. 
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The shock motions induced by the initial perturbations 
are calculated in a finite difference fashion. Suppose that 
the shock location r s h(t n ) and velocity v s h(t n ) are given 
at t n . Then the shock location at the next time £ n +iis 
given by 

r s h(tn+i) = r 3h (t n ) + v 3h (t n ) x At, (1) 

where At = t n+ i — t n is the interval between the two 
successive times. In this study we set At = 10 _4 s, which 
is sufficiently short. The initial time corresponds to n = 

0. If the shock velocity at t n +i is obtained somehow, 
then the procedure is iterated until a designated time is 
reached. 

2.1.4. Step 4- determination of shock velocity - setting up 
Riemann problem - 

How to obtain the new shock velocity v s h(tn+i) at ^n+i 
is the most important part of this model. The idea here 
is similar to the one of the Godonov scheme in numeri- 
cal hydrodynamics: a Riemann problem is set up at the 
beginning of each time step and the evolution during the 
subsequent interval is determined by the solution of the 
Riemann problem. In our model, we consider a certain 
Riemann problem at the new shock location r s h{t n +i) 
and solve it to find the velocity of the forward shock 
that exists always in the solution (see Step 5). 

In order to set up the Riemann problem, we need to 
specify hydrodynamical quantities on both sides of the 
discontinuity located at r s / l (t n+ i). The upstream quan- 
tities are easily obtained from the unperturbed steady 
accretion flow. The downstream quantities are not so 
easy to obtain. In this model we assume that the down- 
stream flow between r s h(t n ) and r a h(t n +i) is steady. We 
solve Eqs. (|A1|) - (|A4|) for steady flows from r s h(t n ) to 
r s h(t n +i) to obtain the hydrodynamical quantities at the 
latter point and use them for the Riemann problem. 

The post-shock flows are not steady. Since they are 
subsonic, it takes the advection time plus sound cross- 
ing time between the shock and the proto-neutron star 
to reach a steady state. In the shock revival phase, how- 
ever, this time scale is longer than the typical time scales 
of variations in the shock radius, which means that the 
steady state is never realized over the entire post-shock 
region. On the other hand, we find in hydrodynami- 
cal simulations that the post-shock flows are approxi- 
mately steady in the vicinity of the shock wave (see Ap- 
pendix [Bj. By virtue of this approximation, we can avoid 
hydrodynamical simulations in the entire region and in- 
stead treat the shock motion alone. 

2.1.5. Step 5: determination of shock velocity - solving 
Riemann Problem - 

Now that the Riemann problem has been constructed, 
we can obtain the new shock velocity, based on the so- 
lution of the Riemann problem. The bottom panels in 
Figure [T] show the schematic pictures of two possible so- 
lutions of the Riemann problem. In the current situa- 
tion, the solution always contain a forward shock wave, 

1. e., the right-going shock wave in Figure [TJ which we 
regard as the displaced shock wave. The left-going wave 
can be either a rarefaction wave (left panel) or a shock 
wave (right panel). Regardless, we adopt the velocity of 
the right-going shock wave as the new shock velocity at 
tn+i, i- e ? v s h(tn+i)- This closes a single time step. We 



go back to Step 3 and repeat the following steps for the 
next time step. The iteration is terminated when the 
designated time is reached. 

2.2. Miscellaneousness 

In this study we fix the neutrino luminosity, mass ac- 
cretion rate and mass of a central object during the time 
evolution for simplicity. It should be stressed, however, 
that the semi- dynamical model can handle the time de- 
pendence of these parameters with no difficulty. It is also 
important to note that the only difference between our 
model and full (ID) hydrodynamical simulations is how 
to obtain the hydrodynamical quantities just behind the 
shock wave. Hence the accuracy of our model depends 
entirely on how good the approximation of local steady 
flows is. As will be demonstrated in $3] and Appendix [Bj 
it looks fine in general. Since the post-shock flows are 
subsonic, they are always affected by the physical condi- 
tions of inner regions in principle. This will be particu- 
larly so if the relaunched shock wave is again stalled. On 
the contrary, if the shock continues to propagate outward 
briskly without stagnation, we expect that our approxi- 
mation will work reasonably well. 

3. RESULTS 

In this section we apply the semi-dynamical method 
to shock revival in the post-bounce supernova cores. We 
first study the characteristics of solutions and observe 
the existence of the critical fluctuations. Then we inves- 
tigate the dependence of the critical fluctuation for shock 
revival on some key parameters, which is our main result 
in this paper. 

3.1. Characteristics of shock evolutions 

We first investigate the evolutions of shock wave, ini- 
tially adding perturbations of different amplitudes to 
the fiducial background model, which is characterized 
by the following conditions: L52 = 5, M S un = 1 and 
Mi n = 1.4M (see section [2.1. ip . Here we define L52 = 
L^/(10 52 erg/s) and M sun = -M/(lM /s) with L v , M 
and M in denoting the neutrino luminosity, mass accre- 
tion rate and mass of a central object, respectively. In 
this fiducial model the standing shock wave is located at 
r s h = 1.5 x 10 7 cm. 

Figure [2] shows three representative evolutions of shock 
radius (left panel) and velocity (right panel). As demon- 
strated in the left panel, the shock wave either settles 
down at a new position (red line) or keep going outward 
(green and blue lines). The green line is in fact a divid- 
ing line of the two cases. These lines correspond to dif- 
ferent initial shock velocities, with the green and black 
lines having larger initial shock velocities. We refer to 
the initial perturbation for the green line as the critical 
shock velocity. The value of the critical shock velocity is 
v sh{crit) — 1-65 x 10 9 cm/s in this case. 

As is evident from the right panel, the out-going shock 
wave is decelerated initially in all cases. The red line, 
which is given half the critical shock velocity initially, 
remain decelerated until the shock wave comes to a halt 
at t ~ 40ms. In the case of the green and blue lines, 
the latter of which is given twice the critical shock veloc- 
ity at the beginning, however, the shock wave begins to 
accelerate at some point of time (t ~ 0.2s for the green 
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Fig. 2. — The evolutions of shock wave for the fiducial background model. The left panel shows the shock radii as a function of time 
whereas the right panel presents the corresponding shock velocities. Three lines with different colors in each panel correspond to different 
shock velocities added initially. See the text for more details. 

line and t ~ 0.01s for the blue line) and keep propagat- 
ing outward up to the end of calculation. It is also clear 
that the shock wave evolves faster for larger initial shock 
velocities among models that do not fizzle out. 

The deceleration is efficient until the shock wave 
reaches r s h ~ 10 8 cm (see Figure [2]). In fact, it seems 
that the shock waves that cannot make this distance 
fail to revive and vice versa. The shock deceleration in 
the early phase of shock revival is consistent with the 
fact that the fiducial background model is linearly stable 
against radial perturbations. Note that the linear sta- 
bility of spherically symmetric shocked accr etion flows 
can be judg ed by the Nakayama's criterion (|Nakavamal 
1 1994 1 19961 ). which states that such flows are linearly 
stable if the post-shock ma tter is decelerated (see also 
lYamasaki fc Yamadal (|2007l )). We confirm that this is 
indeed the case for the fiducial model employed here. Of 
the three models shown in Figure [2j the failed case (red 
line) has relatively small initial perturbations, to which 
linear analysis will be applicable. To the other two cases, 
which are given rather large perturbations, linear analy- 
sis may not be applicable. We think, however, the cause 
of the shock deceleration can be understood in the same 
way at least qualitatively. It is then nice that the semi- 
dynamical model is consistent with linear stability. 

The shock deceleration an d its strength may beun- 
derstood yet another way. iNobuta fc Hanawal (|1994f ) 
claimed that the stability of standing shock wave can 
be understood by the momentum flux balance between 
the pre- and post-shock flows. Of course they are ex- 
actly equal to each other for the initial standing shock 
wave. If the shock wave is slightly dislocated from the 
original position, then the balance is shifted and the 
shock wave will be either pulled back or pushed fur- 
ther. The former occurs if the momentum flux in the 
pre-shock flow is larger than that in the post-shock flow. 
This implies the stability of shock wave. The latter 
case corresponds to instability, on the other hand. Note 
that this stability c riterion is consistent with t hat o f 
Nakayama's (see e.g. iNagakura fc Yamadal (|2008ll2009h : 



iNagakura fc T akahashi (2010)). It is also expected that 
the greater the imbalance of the momentum fluxes is, 
the stronger the deceleration or acceleration will be. 
We apply this criterion to the shock deceleration at 
r s h < 10 8 cm. 

We show in Figure [3] the momentum fluxes just ahead 
of and behind the shock wave if it were dislocated up to 
the specified position by the initial perturbation. The 
three panels correspond to different back groun d mod- 
els (cf. Figure 10 in INobuta fc Hanawal (|1994T >). The 
horizontal axis denotes the radius normalized by the un- 
perturbed shock radius rw ), whereas the vertical axis 
expresses the momentum flux normalized by the value at 
r sh(o)- The red (green) lines give the momentum fluxes 
just behind (ahead of) at the shock wave. The post-shock 
momentum-flux is calculated by extending the steady 
post-shock flow up to the perturbed shock front. On the 
other hand, the pre-shock momentum flux is obtained 
by assuming that the pre-shock flow is unaffected by the 
perturbation. 

The left panel of Figure [3] corresponds to the fiducial 
background model. As is clear from the figure, the exact 
momentum balance is satisfied at the original shock loca- 
tion whereas the pre-shock momentum flux overwhelms 
the post-shock momentum flux as the shock wave is dis- 
located outward. In this case, as already mentioned, the 
shock wave will tend to be pulled back to the original po- 
sition, i.e, the shock wave will experience deceleration. 
We also note that the strength of deceleration differs 
among the background models. The middle and right 
panels of Figure [3] correspond to the cases either with 
a higher neutrino luminosity or with a lower mass ac- 
cretion rate than the fiducial case. In both cases, the 
unperturbed shock waves are located at larger radii than 
in the fiducial model. As is evident in these panels, the 
momentum flux imbalance between the pre- and post- 
shock flows is smaller than in the fiducial model. This 
indicates that the pull back force that decelerates the 
shock wave wanes with the original shock radius. This 
is very important in analyzing the critical fluctuation for 
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Fig. 3. — Momentum fluxes just ahead of and behind the shock wave if it were dislocated to the position specified by the the horizontal 
axis. The shock radius is normalized by the original radius r sh ^ y The vertical axis denotes the momentum flux normalized by the value 
at 7* s h(o)- The red solid line gives the post-shock value whereas the green dotted line presents the pre-shock value. The left panel shows 
the momentum flux imbalance for the fiducial background model whereas the other two panels represent different background models with 
the neutrino luminosity and mass accretion rate displayed in each panel. Both models have larger shock radii initially than the fiducial 
model: r sh ^ = 9.8 x 10 7 cm for the model in the middle panel whereas r sh ^ = 6.4 x 10 7 cm for the model in the right panel. 



shock revival in the next section. 

3.2. The Critical Fluctuation 

In the previous section, we show the existence of the 
critical shock velocity for the fiducial model, i.e., if the 
shock velocity administered initially is larger than this 
value, the shock continues to propagate outward even if 
the neutrino luminosity is smaller than the critical value. 
In this subsection, we perform a larger number of calcu- 
lations for different background models, determining the 
critical shock velocity for each model, and then we ana- 
lyze their property in detail. 

We prepare the unperturbed steady accretion flows as 
follows: we fix the mass of a central object to M^ n = 
1.4M ; in one sequence we vary mass accretion rate 
with a fixed neutrino luminosity (the red line in the left 
panel of Figure |4]) whereas we fix the mass accretion rate 
and change the neutrino luminosity in another series (the 
green line in the same figure). Note that in each sequence 
of background models the location of shock wave (r s h) 
corresponds one-to-one to the free parameter, i.e., either 
the mass accretion rate (for the red line) or the neutrino 
luminosity (for the green line). 

Having in mind that the initial perturbations are gen- 
erated by the hydrodynamical instabilities in reality, we 
refer to the dimensionless perturbation in the post-shock 
pressure, which is defined to be 

/ = (p-Po)M), (2) 
as the fluctuation hereafter. In this expression p and po 
stand for the post-shock pressures for the perturbed and 
unperturbed flows, respectively. The value of the fluc- 
tuation that corresponds to the critical shock velocity is 
called the critical fluctuation. Figure |4] shows the criti- 
cal fluctuations f cr u as a function of the initial radius of 
shock wave (r s h) for the two series of background mod- 
els mentioned above. The critical fluctuation for each 
background model is obtained by many trials of semi- 
dynamical calculations with various initial shock veloci- 
ties. 

As we can see from the left panel of Figure |4j the criti- 
cal fluctuation becomes smaller with the shock radius. It 
is interesting that the two lines are almost identical in the 



range of 10 7 cm < r < 6 x 10 T cm. This is the most impor- 
tant region in discussing shock revival, since it roughly 
corresponds to the typical location of the stagnated shock 
wave. The above results suggest that the critical fluctua- 
tion is mainly determined by the shock location and that 
the neutrino luminosity and mass accretion rate do not 
directly dictate shock revival; instead, they are indirectly 
important in the sense that they determine the initial lo- 
cation of standing shock wave. The reason why large 
radii of standing shock wave are advantageous for shock 
revival can be understood from the momentum-flux im- 
balance. As explained in § I3.ll and Figure [3j the pull 
back force exerted on the shock wave is weakened with 
the radius of standing shock wave and hence relatively 
small fluctuations are sufficient for shock revival at large 
distances. 

It i s interesting to note that Ya masaki fe Yamadal 
(|2QQ5f ) observed almost the identical shock radius for var- 
ious models at their critical points, which suggests that 
the shock location is the important quantity t o determine 
the critical point. iMurphv fc Burrowsl (|2QQ8f ) conducted 
experimental simulations in multi-dimensions and found 
that the larger neutrino luminosity induces the greater 
amplitude of shock oscillations, which is also consistent 
with our interpretation, since the high neutrino luminosi- 
ties result in large shock radii and the pull back force 
acting on the shock wave is then weak. 

We next vary the mass of central object to see the de- 
pendence of our findings on the change. The results are 
shown in the right panel of Figure HI This time we inves- 
tigate the sequence of steady accretion lows, in which the 
mass accretion rate is fixed (M sun = 1) and the neutrino 
luminosity is varied. The result would not be changed 
if the sequence with the fixed neutrino luminosity were 
studied. The red line in the figure represents the criti- 
cal fluctuations for M in = 1.4M (identical to the green 
line in the left panel of Figure 2]) whereas the green line 
denotes the critical fluctuations for M^ n = 1.2M . Al- 
though the difference between the two cases is not so 
large, the critical fluctuation for Mi n = 1.4M is sys- 
tematically larger by several percentage points than that 
for Mi n = 1.2M for the same initial shock radius. This 
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Fig. 4. — The critical fluctuations as a function of the initial shock radius. The horizontal axis denotes the original shock radii in the 
unperturbed steady accretion flows, whereas the vertical axis represents the normalized critical fluctuations, which is defined by Eq. ([2]). 
In the left panel the red line represents the sequence with a fixed neutrino luminosity (L52 = 5) whereas the green line corresponds to the 
series with a fixed mass accretion rate (M S un = 1)- The mass of the central object is set to Mi n = IAMq. In the right panel, on the other 
hand, we compare the results for different masses of the central object. The red line represents the model with Mi n = IAMq whereas the 
green line corresponds to the model with M^ n = 1.2Mq. This time the mass accretion rate is fixed (M S un = 1) and the neutrino luminosity 
is varied. 



is mainly attributed to the fact that gravity is stronger 
for the heavier central object, which then leads to the 
greater ram pressure. The result clearly demonstrates 
that larger masses of the central core are negative for 
shock revival, which is consis tent with the analy s is bas ed 
on the critical luminosity by iKeshet fc Balbergl (|2Q12[ ) . 

We find that the results given in Figure [4] can be ap- 
proximated by 

fc^t ~ 0.8 x (J^L.) x {1 - (-^-)}. (3) 
Jcrit 4.4M© ; 1 V 10 8 cm y/ V J 

This simple analytic formula will be useful in analyz- 
ing the onset of explosion in full dynamical simulations, 
since the fluctuation / can be easily estimated at each 
time step. Note that po is obtained from the values of hy- 
drodynamical quantities just ahead of the stalled shock 
wave by using the Rankine-Hugoniot relation for v s h = 0. 
We expect that if the fluctuation / at a certain time step 
in a simulation does not exceed f cr u estimated this way, 
the shock wave does nothing but oscillations around the 
average position. 

The critical fluctuation given by Eq. (|3)) might be a 
conservative estimate, since we ignore the time evolution 
of mass accretion rate in this analysis. Since the mass 
accretion rate becomes smaller with time in reality, the 
ram pressure at shock front, which is one of the main ob- 
stacles for shock propagation, should be smaller, which 
would likely help the shock wave move outward. In fact, 
we confirm in another series of semi-dynamical calcula- 
tions that decrease in the mass accretion rate can induce 
shock revival (see Appendix [CJ . 

4. SUMMARY AND DISCUSSIONS 

In this paper, we have developed a new approach to 
shock revival in the neutrino heating mechanism in core- 
collapse supernova theory. The semi-dynamical model 
takes into account the temporal evolution of shock wave 
induced by fluctuations in the hydrodynamical instabili- 



ties. The equation of motion of the shock wave is mod- 
eled in a finite difference manner, setting up a Riemann 
problem at each time step. In so doing the post-shock 
flows are approximated by locally steady flows, which 
are indeed observed in hydrodynamical simulations. The 
method is very efficient in conducting a large number of 
model calculations of the long-term shock evolution. 

As a result, we have demonstrated that there is a crit- 
ical fluctuation for a given steady accretion flow, above 
which the shock can sustain outward propagations, lead- 
ing hopefully to a successful supernova explosion. Vary- 
ing the neutrino luminosity or mass accretion rate for 
different massed of a central object, we have obtained 
a simple fitting formula for the critical fluctuation as a 
function of the initial shock radius and mass of central 
object (see Eq. fl3J). According to our results, the crit- 
ical fluctuation decreases with an increase in the shock 
radius and/or decrease in the mass of central object. In 
particular, we have found that the initial shock radius is 
the key parameter for shock revival. 

It should be noted that in our models the neutrino 
luminosities are all sub-critical and there is a steady 
shocked accretion flow for a given pair of neutrino lumi- 
nosity and mass accretion rate. In fact, even if the shock 
wave is stagnated around r s h ~ 10 7 cm, in which case 
the neutrino luminosity is much smaller than the criti- 
cal value, a large enough fluctuation can put the shock 
wave into a sustained outward propagation. We have 
hence concluded that although the neutrino luminosity 
and mass accretion rate are important in determining 
the initial location of stagnated shock wave, it is the dy- 
namical fluctuations that have direct leverage on shock 
revival. 

Based on our findings in this paper, we now discuss 
possible implications for the neutrino heating mecha- 
nism in multi-dimensional hydrodynamics. According 
to the recent results of 3D core-collapse simulations by 
iTakiwaki etHI (|2012t ). the maximum residency time of 



8 



Nagakura et al. 



1e+08 



CO 



1e+07 



I I I 


1 1 

/ • ■ ■ 
/ 1 • i 
/ 1 ■ ; 
/ 1 * j 




/ ; 




1 . 






L52 =j - 




: £52 =_3 - 


p-i-r 


11 1 





1e+08 



^ 1e + 07 





* 1 J 
c ' 1 / 
: • f 1 

ill 
1 / 

1 / ~ 




* * / 




*** '* 




M SU n = l : 




M aun = 0.5 : 




M sun = 0.1 
■ 1 



-2.5 -2 -1.5 -1 

M (M /s) 



-0.5 



1e+52 



L(erg/s) 



1e+53 



Fig. 5. — Radii of standing shock waves in unperturbed steady accretion flows. In the left panel, we vary the mass accretion rate for 
fixed neutrino luminosities whereas in the right panel, we fix the mass accretion rate and change the neutrino luminosity for each line. 



accreting matter in the gain region is longer in 3D than 
in 2D owing to the extra degree of freedom in motion. 
This leads to longer neutrino heating, which then may 
result in larger average radii of stagnated shock wave in 
3D than in 2D. According to our findings in this paper, 
the larger shock radius means a smaller critical fluctu- 
ation for shock revival. In this sense at least 3D mod- 
els are more favorable for successful explosions. How- 
ever, the dynamical fluctuations generated in SASI or 
neutrino-driven convections may be smaller in 3D than 
in 2D, since the free energy of turbulence could be dis- 
trib uted over a larger number of oscillation modes in 
3D (|Nordhaus et aLlhoiOl : iBurrows et al.ll2012l ). In fact, 
the sloshing modes, which are always observed markedly 
in 2D axisymmetric simulations, are not so remarkable 
in 3D with non- axisymmetric modes being dominant. 
These two factors (the initial shock location and the 
amplitude of fluctuations in hydrodynamical instabili- 
ties) compete with each other and make it difficult to 
unravel the net effect of 3D hydrodynamics, the fact 
that may lead to the current contr o versy between dif- 
ferent groups (iNordhaus et all [201(1 lHanke et al"1l20lH 
iTakiwaki et aLll2Q12[ ). 

Although it is beyond the scope of this paper to set- 
tle the dispute on the 3D effect, we can add some more 
considerations that may be useful. Figure [5] shows the 
locations of standing shock wave in unperturbed steady 
accretion flows as a function of mass accretion rate (left 
panel) or neutrino luminosity (right panel) . In the former 
the neutrino luminosity is fixed whereas the mass accre- 
tion rate is fixed in the latter. As shown clearly for both 
cases in this figure, the location of standing shock varies 
very rapidly as the mass accretion rate (left panel) or 
neutrino luminosity (right panel) exceeds a certain value. 
For instance, the blue line in the left panel, for which the 
neutrino luminosity is fixed at L52 = 5, shows that r s h is 



almost constant as long as |M snn | > 1 whereas it grows 
rapidly for |M snn | < 0.5. A similar trend is evident in the 
right panel if one replaces the mass accretion rate with 
the neutrino luminosity. This may imply then that if the 
average shock radius is already in this rapidly changing 
regime, the shock radius of the two competing factors 
will be more important and 3D hydrodynamics may be 
more advantageous for shock revival. If the shock wave 
is stagnated at a rather small radius, on the other hand, 
the larger fluctuations in 2D will be more important to 
provide more favorable conditions for shock revival. 

In spite of its simplicity the semi-dynamical approach 
we have employed in this paper can capture the essential 
features of the shock dynamics such as linear stability. 
All the results we have obtained are of qualitative nature 
and the critical fluctuations will be somewhat changed if 
more realistic treatments of hydrodynamics and micro- 
physics are incorporated. We believe, however, that the 
existence of the critical fluctuation and its qualitative de- 
pendence on the shock radius and mass of central object 
will not be changed. It is true that there are some limi- 
tations to the semi- dynamical approach in this paper. In 
addition to the fact that it is ID, some features in shock 
dynamics, such as the evolution after re-stagnation, cer- 
tainly cannot be described, since the locally steady ap- 
proximation will not be valid. More detailed compari- 
son with multi-D hydrodynamical simulations will reveal 
both the merit and demerit of the approach, which will 
be the issue of our forth coming paper. 
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APPENDIX 

APPENDIX A: CONSTRUCTION OF STEADY SHOCKED ACCRETION FLOWS 

In this section, we present details in the construction of unperturbed, spherically symmetric, steady accretion flows 
in Step 1 of §2.1.11 We solve the following ordinary differential equations from the shock front (r s h) to the neutrino 
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Fig. 6. — Comparison of the radial profiles of density (left panel) and pressure (right panel) between the hydrodynamical simulation 
model and the steady accretion flow. See the text for details. 

sphere for electron- type neutrinos (r Ue ): 

Anr 2 pv M, (Al) 

dv 1 dp GMin n , ao v 

v— + -^- + — 5— =0, (A2) 
dr p dr r z 

d e p dp Q E , aq x 
dr p p A dr p 

= Q N . (A4) 

dr 

In the above expressions, p, p, T, e and Y e denote the mass density, pressure, temperature, energy density and electron 
fraction, respectively. Other symbols, r, v, G, M, Mj n , stand for the radius, fluid velocity, gravitational constant, 
mass accretion rate and mass of a central object, respectively. Here Newtonian gravity is assumed and the self-gravity 
of accreting matter is neglected. Interactions between neut rinos and matter are encapsulated in Qe and Qn, the 
expressio ns of which are ad opted from Eqs. (16) and (17) in lOhnishi et al.l (j2QQ6h . We employ the equation of state 
(EOS) by Sh en et aTl (|201l[ ). which is based on the relativistic mean field theory. 

At the shock surface, we impose the Rankine-Hugoniot relation for the standing shock wave. The pre-shock flow is 
assumed to be a free fall with the entropy per baryon of s = 3ks and electron fraction of Y e = 0.5. Here ks is the 
Boltzmann constant. At the inner boundary or the neutrino sphere (rv), on the other hand, the following relation is 
assumed: 

Lu = ^T^rl, (A5) 

where L V1 a and T v denote the neutrino luminosity, Stefan-Boltzmann constant and neutrino temperature, respectively. 
We further assume that the neutrino luminosities for electron- type neutrinos (L Ue ) and anti- neutrinos (L^) are identical 
(and denoted as L v ). We ignore other types of neutrinos, which play only a minor role in matte heating. For all the 
models investigated in this paper, the neutrino temperatures are fixed to T Ue = 4MeV and T„ e = 5MeV, the typical 
values in the post-bounce phase. Note that given L v , the neutrino spheres (r Ue and rv e ) are uniquely determined in 
this prescription. 

Once the shock radius (r s h), mass of the central object (M in ), neutrino luminosity [L v ) and mass accretion rate 
(M) are specified, a steady post-shock accretion flow is uniquely obtained as a solution of Eqs. (|Al|) - (jA4|) . Imposing 
further that the density at the inner boundary be p = 10 11 g/cm 3 , we determine the shock radius r s ^. The results 
given in Figure [5] are also obtained this way. 

APPENDIX B: VALIDITY OF THE LOCAL APPROXIMATION 

The semi-dynamical model assumes that the post-shock flows are locally steady near the shock wave. As mentioned in 
the text, the reliability of the approach depends cr itically on the validity of this assumption. Here we show a comparison 
with a result of hydrodynamical simulations bv lYamamoto et al.l (|2Q12f ). which may support our assumption. The 
hydrodynamical models employ the so-called light-bulb approximation. We pick up a model of their computations, 
which has following parameters: L52 = 5, M sun = 0.53 and Mi n = 1.4M . We adopt the same values for these 
parameters also in the semi-dynamical model in this comparison. 
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Fig. 7. — The time evolutions of shock radius (left) and velocity (right) after reducing the mass accretion rate by 1%. In the upper 
panels, the background models have the same neutrino luminosity (L52 = 5) but different mass accretion rates, whereas the background 
models for the lower pannels have the identical mass accretion rate (M SU n = 1) but different neutrino luminosities. 

We take a snap shot from the hydrodynamical data and compare it with the corresponding steady post-shock 
configuration. For that purpose we need to extract the shock velocity in the hydrodynamical simulation. We first 
identify the pre- and post- shock states and then calculate the shock velocity from the momentum flux balance. 
Although the numerical diffusion in dynamical simulations makes it difficult to identify the pre- and post-shock states 
very accurately, we decide to adopt the hydrodynamical quantities at the points, where the density gradient starts to 
change rapidly. Although this is a rather crude treatment, it is sufficient for the present purpose. 

Figure [H shows the results of comparison. The left panel in this figure displays the density distributions whereas the 
right panel gives the pressure distributions. The red lines labeled as SD represent the steady flow whereas the green line 
with the label FD gives the results of the dynamical simulation. The shock velocity calculated from the momentum 
balance is ~ 1.3 x 10 9 cm/s in this case, which exceeds the critical shock velocity obtained by the semi-dynamical 
approach and, in fact, lead to shock revival in the dynamical simulation. As shown in this figure, the density and 
pressure distributions near the shock front are almost identical with those of the steady accretion flow. We hence 
confirm at least in this case that the approximation of local steadiness is good. 

APPENDIX C: SHOCK EVOLUTIONS INDUCED BY THE REDUCTION OF MASS ACCRETION RATE 

Although we fix the mass accretion rate for simplicity in this paper, it is a function of time in reality and is an 
important factor to induce motions of stagnated shock waves. In this appendix, we consider the shock evolution 
induced by the reduction of mass accretion rate imposed initially by hand instead of taking fully into account the 
time-dependent mass accretion rate. 

After setting up a steady accretion flow (Step 1), we artificially reduce the mass accretion rate in the pre-shock region 
by a certain fraction. This naturally leads to a Riemann problem at the standing shock wave and induces its motion 
without perturbations to the shock velocity. The subsequent steps in the semi-dynamical method are completely 
the same as before (Section [2]). Note that the change of M is given only at the beginning and it remains constant 
afterwards. It is also assumed that the neutrino luminosity and mass of a central object {Mi n = 1AM@) are constant 
in this study. 
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Fig. 8. — The time evolutions of shock radius (left) and velocity (right) after reducing the mass accretion rate by different fractions. 
The neutrino luminosities and mass accretion rates in the background models are given in the panels. The red, green, blue and pink lines 
correspond to the initial reductions by 1, 10, 50 and 90% of mass accretion rate, respectively. 

Figure shows the evolution of shock radius (left panels) and velocity (right panels) for the models, in which M is 
reduced by 1% initially. The vertical axis of left panels is the shock radius normalized by the initial value, whereas the 
horizontal axis denotes the time after the onset of calculation. The upper panels present the results for the background 
models with a fixed neutrino luminosity of [L v = 5 x 10 52 erg/s) but different mass accretion rate, whereas the lower 
panels display the results for the background models with a fixed mass accretion rate [M sun = 1) but with different 
neutrino luminosities. 

In all models, the reduction of mass accretion rate in the pre-shock flow leads to the outward motions of the shock 
wave just as expected from the initial Riemann problems. Intuitively speaking, the smaller mass accretion rate implies 
weaker ram pressures. As the neutrino luminosity increases or the mass accretion rate decreases in the background 
model and, as a consequence, the initial shock radius becomes larger, the shock reaches greater distances later. This 
trend is similar to the case considered in the main body, in which velocities are added to the standing shock wave 
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initially. This is understood from the pull back forces exerted on the shock wave just in the same way. 

Next we change the initial reduction rate in the mass accretion rate of the pre-shock flows. Figure [8] shows the 
evolution of shock waves for 3 different background models. The upper panels correspond to the fiducial background 
model, whereas the middle (lower) panels show the results for the background model with M sun = 1 and L52 = 10 
(M SU n — 0.2 and L52 = 5). The left panels present the evolutions of shock radius and the right panels display 
the temporal changes of shock velocity. The red, green, blue and pink lines correspond to 1%, 10%, 50% and 90% 
reductions in the mass accretion rate, respectively. As is expected again, larger changes in the mass accretion rate 
pushes the shock wave outward more strongly (see Figure [8j). In fact, we can consider the critical reduction rate in 
the mass accretion rate, below which the shock wave continues to move outward without re-stagnation. 

The critical reduction of mass accretion rate may be useful when one judges whether a hit of an interface between 
different layers, e.g., Si/O interface, on the stagnated shock wave will lead to shock revival. It is also clear that a secure 
decline of mass accretion rate, which is ignored in this paper, will be helpful for the shock propagation induced by the 
fluctuations in hydrodynamical instabilities. In this sense the critical fluctuations given by Eq. (j3j) is a conservative 
estimate. 
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